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We calculate the geometric phase for an open system (spin-boson model) which interacts with 
an environment (ohmic or nonohmic) at arbitrary temperature. However there have been many 
assumptions about the time scale at which the geometric phase can be measured, there has been 
no reported observation of geometric phases for mixed states under nonunitary evolution yet. We 
study not only how they are corrected by the presence of the different type of environments but 
estimate the corresponding times at which decoherence becomes effective as well. These estimations 
should be taken into account when planning experimental setups to study the geometric phase in 
the nonunitary regime, particularly important for the application of fault-tolerant quantum compu- 
tation. 
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I. INTRODUCTION 

All real world quantum systems interact with theirs 
surrounding environment to a greater or lesser extent. 
While closed quantum systems are bound to have an 
unitary evolution in which the system's purity is pre- 
served and the superposition principle can be applied, 
open quantum systems show a different scenario. As they 
are in interaction with an environment (defined as any de- 
grees of freedom coupled to the system which can entan- 
gle its states), a degradation of pure states into mixtures 
takes place. These mixtures states will often turn out to 
be diagonal in the set of "pointer states" which are se- 
lected by the crucial help of the interaction Hamiltonian 
{Hint)- They are stable subjected to the action of Hint, 
i.e., the interaction between the system and the environ- 
ment will leave them unperturbed. That's exactly what 
makes them a "preferred" basis. No matter how weak the 
coupling that prevents the system from being isolated, 
the evolution of an open quantum system is eventually 
plagued by nonunitary features like decoherence and dis- 
sipation. Decoherence, in particular, is a quantum effect 
whereby the system loses its ability to exhibit coherent 
behaviour. Nowadays, decoherence stands as a serious 
obstacle in quantum information processing. 

Since the work of Berry the notion of geometric 
phases (GPs) was shown to have important consequences 
for quantum systems. Berry demonstrated that quantum 
systems could acquire phases that are geometric in na- 
ture. He showed that, besides the usual dynamical phase, 
an additional phase that was related to the geometry of 
the space state was generated during an adiabatic evolu- 
tion. The original idea, framed within the context of adi- 
abatic and cyclic evolutions of isolated systems, has been 
generalized in various aspects. While many of this propo- 
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sitions have been centered around pure states, the need to 
address the issue of GPs for mixed states rapidly gained 
prominence fuelled by the promise of realizing quantum 
logic gates under realistic physical conditions. In that di- 
rection, in 01, it was introduced an alternative definition 
of GPs for nondegenerate density operators based upon 
quantum interferometry. In Q, a kinematic description 
of the mixed state GP was given and its definition ex- 
tended to degenerate density operators. Recently, in j^, 
GPs for an open quantum system were studied and, fi- 
nally, it was shown that the above apparently different 
approaches were related in a unifying framework. The 
effect on GPs of differents types of decoherence sources, 
such as dephasing and spontaneous decay, has been an- 
alyzed 0. Likewise, it has been shown how to generate 
a geometric phase through modifications solely on the 
reservoir that interacts with a small subsystem 01 . 

GPs are useful in the context of quantum computation 
as a tool to achieve fault tolerance. However, practi- 
cal implementations of quantum computing are always 
done in the presence of decoherence. Thus, a proper gen- 
eralization of the geometric phase for unitary evolution 
to that for nonunitary evolution is central in the eval- 
uation of the robustness of geometric quantum compu- 
tation. This generalization to nonunitary evolution has 
been done in Q , where a functional representation of GP 
was proposed, after removing the dynamical phase from 
the total phase acquired by the system under a gauge 
transformation. 

The GP for a mixed state under nonunitary evolution 
is then defined as 

$ = arg{;^ ^Ek (O)efe (r) (0) | (r)) 

k 

X g-Jo-dt(*.|i^|*.>}^ (1) 

where ek{t) are the eigenvalues and j^'fc) the eigcnstates 
of the reduced density matrix (obtained after tracing 
over the reservoir degrees of freedom) . In the last defini- 
tion, T denotes a time after the total system completes 
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a cyclic evolution when it is isolated from the environ- 
ment. Taking the effect of the environment into account, 
the system no longer undergoes a cyclic evolution. How- 
ever, we will consider a quasicyclic path V : t e [0,t] with 
T = 27r/i7 (n the system's frequency) It is worth not- 
ing that the phase in Eq.(^ is manifestly gauge invari- 
ant, since it only depends on the path in the state space, 
and that this expression, even though is defined for non 
degenerate mixed states, corresponds to the unitary ge- 
ometric phase in the case that the state is pure (closed 
system) 

It is expected that GPs can be only observed in exper- 
iments carried out in a time scale slow enough to ignore 
nonadiabatic corrections, but rapid enough to avoid the 
destruction of the interference pattern by decoherence 
So far, there has been no experimental observation of 
GPs for mixed states under nonunitary evolutions. The 
purpose of this letter is to study how GPs are affected 
by decoherence. Not only shall we analyze the effect of 
the environment on the GPs and their robustness against 
decoherence but also, under which conditions GPs can be 
measured. With these motivations, we shall introduce a 
"simplified" spin boson model and calculate the correc- 
tions to the unitary geometric phase for different physical 
environments. In the end, we shall estimate the decoher- 
ence times after which GPs can be no longer measured 
because they, literally, disappear. 

The paper is organized as follow. In Section II we 
introduce the spin-boson model and we show the exact 
master equation for the reduced density matrix. We ex- 
plicetly evaluate, in Section III, the corrections to the 
GPs induced by different types of environment and Sec- 
tion IV contains the calculation of the decoherence times. 
Finally, in Section V, we make our final remarks and com- 
ments. 



II. THE SPIN-BOSON MODEL 

Hereafter, we shall study an exact model for a two- 
state quantum system [or quantum bit (qubit)] coupled 
to a thermal bath of harmonic oscillators, where deco- 
herence is the only effect on the system-particle. This 
is a particular case of the spin-boson model of Ref.[lol| 
(where the tunnelling bare matrix element is A = 0) and 
has been used by many authors to model decoherence 
in quantum computers pH and, in particular, it is ex- 
tremely relevant to the proposal for observing GPs in a 
superconducting nanocircuit We will concentrate 

on ohmic and supraohmic environments (at high, low, 
and zero temperature) coupled to the system spin. The 
Hamiltonian that describes the complete evolution of the 
two-state system interacting with the external environ- 
ment is: 

HsB = ^fl^CTz + ^O-j ^ KiQn + ^ f^nO-ian, (2) 
n n 



where the environment is described as a set of harmonic 
oscillators with a linear coupling in the oscillator coordi- 
nate. The interaction between the two-state system and 
the environment is entirely represented by a Hamiltonian 
in which the coupling is only through CTr. A coupling of 
this form indicates that the state of the environment will 
be sensitive to the values of az, which means that the 
environment can "observe" the value of cFz (i.e, where 
the system is | |> or | \,>). The reason for this type 
of coupling is that the effects of a coupling proportional 
to <7x and/or ay can be completely taken into account 
by the renormalization of the natural frequency of the 
system. However, in this particular case, [<T^,i?int] = 
and the corresponding master equation is much simpli- 
fied, with no frequency renormalization and dissipation 
effects. If we assume that (i) the system and the environ- 
ment arc initially uncorrelated and (ii) the environment 
is initially in thermal equilibrium at temperature T (or 
the vacuum for zero temperature), the master equation 
for the reduced density matrix is 

p, = -ir2[cr^,pi.] - 'D{t)[(T^, [<Jz,pA]- (3) 

In other words, the model describes a purely decoher- 
ing mechanism, solely containing the diffusion term 2?(t), 
where no energy exchange between system and bath is 
present. The diffusion coefficient is given by 

V{t)= [ ds [ dujJ{uj)cos{ujs)coth{^^), (4) 
Jo Jo 2 

where (3 = l/ksT (ks is Boltzman constant) and J{uj) 
is the spectral density of the environment defined by the 
expression J(w) = J^n^n^i'-^ ~ ^n)/2'mnUJn- One as- 
sumption we shall make throughout this paper is that 
J{uj) is a reasonably smooth function of w, and that is of 
the form up to some frequency A that may be large 
compared to Q. In particular, the case with n = 1 is the 
"ohmic" case and the one with n > 1 is the "supraohmic" 
case. It is easy to see in Eq.Q that the system's dynam- 
ics depends crucially both on the external temperature 
and on the type of environment. 

III. GEOMETRIC PHASE FOR THE 
COMPOSITE SYSTEM 

We shall compute the geometric phase for the spin- 
boson model. It is easy to check that Proi(t) = 
g-if2t-r(t)^^^^ j-g-j |g solution of the master equation 

(Eq.®) with T{t) = dt'V{t') and Proi(O) a constant 
determined by initial conditions. Here we are imply- 
ing that the diagonal terms do not evolve in time, i.e. 
Prii(i) — Prii(O), for i = 0, 1 (the spin state basis), assum- 
ing the population dynamics of the system is essentially 
frozen on the time scales of interest [ll| . The decoherence 
factor is then given by 

f°° , , , , f Bnw\ (l-cos(wi)) , , 
T{t) ^4 J duj J{uj) cothf -^^ " . (5) 
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When a specific choice turns out to be useful, we shall 
assume the following functional form for the spectral den- 
sity J{uj) = 7o/4a>"A"-ie"'^/^ 113, where 70 is the dis- 
sipative constant (in suitable units) and A is the cutoff 
frequency. 

As the factor T is real, unitarity requires proi (0) = 
/Orio(O) since p-c^„ = Proi*- For a spin 1/2 system, the 
state space consists of all points on and inside the Bloch 
sphere. Then, we assume that the system is initially in 
the Bloch state 



= cos(0o/2)|e) + sin(eo/2)|5). 

The constants are determined so that Pr(0) = 
|«'(0))(«'(0)|. Then, for times t > 0, the reduced den- 
sity matrix is 



Pr(t) = 



cos(0o/2)2 
l/2sin(6'o)e-'"*-r(t) 



l/2sin(6'o)e*f"-r(*) 
sin(0o/2)' 



The eigenvalues of the above reduced density matrix 
are easily calculated: 

e±{t) = i ± iVcos(0o)2 +exp(-2r(t))sin(0o)2. (6) 

In order to estimate the geometric phase, wc only need 
the eigenvector since £-(0) = 0, and, hence, the 

only contribution to the phase comes from that eigenvec- 
tor and its corresponding eigenvalue (see Eq.(^). Then, 
we write it here as: 



-iflt 



sm{9t/2)\e)+cos{9t/2)\g), (7) 



with tan(6't/2) = exp(-r(i)) cot(6'o/2). It is easy to 
check that, when F = 0, we reobtain the results for the 
unitary case. Hence, once we have the eigenvalues and 
eigenvectors, we can calculate the factor (^fclvPfc) of Eq. 
(^. Performing the time derivatives, we get (^'fcl^'fc) = 
-~iilsm{9t/2)^ , therefore, the geometric phase is 



$ = arg v/e+(T)£+(0)(vI'+(r)|vI/+(0))e*"/cr'i*-n(^)= 



This result holds for any density matrix has these eigen- 
values and eigenvectors, independently of the exact ex- 
pression of r, as long as e_ (0) = 0. By inserting Eqs. ® 
and (0 into the above expression, the geometric phase 
related to a quasicyclic path V : t e [0,t] (with r = 27r/r2 
[3) is (assuming cos(6'o/2) > 0): 



<P = n dt sin i — 



(8) 



In order to evaluate this integral, we will perform a se- 
rial expansion in terms of powers of the dissipative con- 
stant 7o, and, consequently, the unitary phase is cor- 
rected by the presence of the environment as: 



$^ + (5$ «7r(l -cos(6lo)) 



7o^sin(6'o)^cos(6'o) 



dt 



am 



where we can see that the first term of the expansion 
is the solution we would have obtained if the evolution 
would have been unitary, i.c — 7r(l — cos(6'o)). If, 
for example, we assume an ohmic (n = 1) environment 
at high temperature, then the spectral density takes the 
particular form J{lu) = 7o/4 Loer^l^. In that case, 
hw << 2kBT, wc can approximate coth(f3huj/2) in Eq. 
Q by 2kBT/ {hui) and the decoherencc coefficient F be- 
comes F = {jQTTkBT)t/h. Then, the geometric phase is 
corrected by: 

This correction to the unitary phase is proportional to 
^qUbT /Ml. Then, it is bigger for hotter environments 
and can not be neglected. It is worth mentioning that, 
formally, our solution is similar to the one proposed by 
authors in where they made no a priori assumption 
about the dynamics of the system. On the contrary, 
our solution is obtained from scratch and can be applied 
to different environments at any temperature. What's 
more, this correction is more realistic since they assumed 
a very much simplified environment. 

If we were to assume the same ohmic environment but 
at zero temperature, then the phase would be corrected 
in a significantly different way, as one would expect since 
there is one time scale lost. In that case, the factor 
coth(/3fiti;/2), in the definition of V{t), (Eq. gjl) can be 
approximated by 1 and the correction obtained is: 



5$ 



7r7o 



l + log(^))sin((?„)2cos(0o)- (10) 



It is worthy noting that this correction comes from the 
zero point fluctuations of the environment. 

Beyond the commonly assumed ohmic spectrum for 
the bath, generic nonohmic environments can be studied 
with this model. Owing to the ultrashort time bath cor- 
relations, nontrivial short-time system dynamics enters. 
What's most appealing about this case is the fact that 
the electromagnetic field can be modeled by a supraohmic 
environment which results very useful for quantum optics 
and trying to measure the geometric phases. In the fol- 
lowing, wc shall evaluate the correction of the geometric 
phases in the presence of this kind of environments. 

For high temperature, we can still make the same ap- 
proximation than in the ohmic case, and then, the es- 
timation of the F factor is straightforward. Using the 
spectral density with n = 3, the correction to the uni- 
tary phase is: 



5$ = 7r7o 



2kBT 

hA 



sin(6'o)^ cos(0o). 



(11) 



We can also estimate the correction to the unitary phase 
in the presence of an supraohmic environment but at zero 
temperature as: 



70=0 



= 7r7o sin(6'o)^ cos(6'o). 



(12) 
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As A is the maximum available frequency in the environ- 
ment, a valid assumption is that f2 < A. In that case, 
we can see that the parenthesis in Eq. 1)10(1 is of order one 
and, then, both corrections at zero temperature fEa. 1)10(1 
and Ea. lll2() ). for the ohmic and nonohmic case, are sim- 
ilar. The same occurs for high temperate in Eq.® and 
Eg. ((11(1 . These results arc remarkably interesting and en- 
hance the robustness of the model (and of the geometric 
phase). 



IV. DECOHERENCE TIMES FOR THE 
COMPOSITE SYSTEM 

In this section we shall estimate the decoherence time, 
i.e. the time-scale when the nondiagonal terms of the 
reduced density matrix are suppressed. Therefore, it is 
a very important time scale to take into account when 
planning an experiment to measure the geometric phases. 

Let's take for example an interference experiment. The 
experiment starts by the preparation of two wave packets 
in a coherent superposition, assuming each of the par- 
ticles follows a well defined classical path (Ci and C2, 
respectively), as ^{t = 0) = ((^i(x) + f2{x)) Xoiv), 
where Xoiv) represents the initial quantum state of the 
environment (whose set of coordinates is denoted by y). 
Due to the interaction between the system and the en- 
vironment, the total wave function at a later time t is 
*(*) = 'Pi{x,t) (g) Xi{y,t) +<f2{x,t) (g) X2{y,t). It is easy 
to note that the states Lpi and 1^2 became entangled with 
two different states of the environment. Therefore, the 
density matrix of the system at time t (for example when 
interference pattern is examined) is, 

Pr{x,x') = Lp\{x)lPi{x') + ipl{x)Lp2{x') 

+ \V*l{x)if2{x) + if*2{x)Lpi(x')]j'{t). 
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FIG. 1: Comparison between the T{t) coefficients for the 
ohmic case (n = 1) and the supraohmic case (n = 3) for 
difi'erent environments for high (T = 1000), low (T = 1.55) 
and zero temperature {h = 1 = ks)- A — WOO, and 70 = 0.3 
(top) and 70 = 0.03 (bottom). Time is measured in units of 
Q. It is easy to note that the low temperature limit is not 
zero temperature. 



The last two terms in the above expression represent 
the quantum interferences. The overlap factor T{t) 
encodes the information about the statistical nature of 
noise since it is obtained by tracing over the degrees 
of freedom of the bath. Hence, it produces a decaying 
term that tends to eliminate the interference pattern. 
Therefore, noise makes J^{t) less than one, and the goal 
is to quantify how it slightly destroys the particle in- 
terference pattern. In the case of the exact system wc 
are studying here, T{t) = e^^'*^ When F is big, the 
off-diagonal terms (coherences) elements will vanish in a 
short time scale, and hence, the open system will not ac- 
quire GP. Then, in order to know the time-scale at which 
the quantum interferences arc suppressed, wc have to es- 
timate the "decoherence" time tjj for each of the four 
cases studied above, by setting T{tD) ~ 1- 

For the ohmic environment, the decoherence factor is 
F^"^ = {kgTTT^Q/fijt in the limit of high temperature. 
Then, the decoherence time scale for this case is, in prin- 
ciple, extremely short, estimated as t]^ = h/ (kBTTT^o) 



(see Fignj. All other cases depend explicitly upon the 
cutoff. For the ohmic environment at zero temperature, 
the decoherence factor is F^° = 7o/21og(l + h?t^). For 
times h.t > 1, the decoherence time scales as t^^ = 
Sooner or later, decoherence is always present 
in this environment as well. As it is clearly shown in 
FigH decoherence is delayed as 70 decreases (bottom). 
Things are slightly different for the supraohmic environ- 
ment. In the high temperature limit, the decoherence 
factor is F^'^ = (2kBT-/o)At'^ /h for At < 1. Then, 
the estimation of the decoherence time is straightfor- 
ward t|f = l/A^M/(2fcsr7o)- However, if At > 1, 
F'^-^ = 2kBT-jo/{hA), constant in time. Then, decoher- 
ence shall occur if and only if 2fcsT7o ^ HA, and that 
shall happen in a time t < 1/A as it is shown on top of 
Figffl At the bottom of Fig^ 70 is smaller and then, 
kBT-fo < HA, which means that F'^-^ will never be of or- 
der one and. therefore, decoherence shall not be effective. 
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Finally, in the supraohmic case at zero temperature, the 
decoherence factor is = JoAH'^/{1 + A'^t^Y . We can 
see in Fig^that decoherence never occurs for this case. 
This is so because F'^'' = 70 for At > 1. The decoher- 
ence factor will be a constant value of e~^° for all times. 
Since 70 < 1, this factor will never be of order one. In 
conclusion, as high-T environments are so efficient pro- 
ducing decoherence, there is a strong condition over the 
dissipative constant 70 in order to get to > t [l5l |. i.e. 
7o < hA/k^T, which actually is unfeasible. Bigger 70 
implies shorter decoherence times. However, for zero- 
T baths, the condition on 70 relaxes to 70 < 1 making 
tu > T and allowing experimental observation. 

V. FINAL REMARKS 

When measuring a GP, the dynamical phase has to be 
eliminated. It can be either canceled, for example, using 
spin-echo technique for spins in magnetic fields or 
one can parallel transport the state vector in order to en- 
sure that the dynamical phase is zero at all times. Even 
though the state does not acquire a phase locally, it can 
acquire a phase globally after completing a cyclic evolu- 
tion. This global phase is equal to the geometric phase. 
GPs are useful in the context of quantum computation as 
a tool to achieve fault tolerance. However, practical im- 
plementations of quantum computing are always done in 
the presence of decoherence, a nonunitary effect of open 
systems. However, there have been many assumptions 
about the time scale at which the geometric phase can 



be measured, there has been no reported observation yet 
for mixed states under nonunitary evolutions. 

It has been argued that the observation of GPs should 
be done in times long enough to obey the adiabatic ap- 
proximation but short enough to prevent decoherence 
from deleting all phase information. In this article, not 
only have we shown an exact model where the correction 
to the GPs can be evaluated for different types of envi- 
ronment at any temperature but also, estimations of the 
corresponding times at which the interference pattern is 
mostly reduced and decoherence becomes effective. This 
model allows us to evaluate an exact master equation 
for the reduced density matrix and, as a consequence, to 
elaborate a full analysis on the effect of decoherence on 
the geometric phase of the composite system. The goal is 
the analysis of the effect produced by decoherence, which 
should be essential to design experimental setups in order 
to observe geometric phases using, for example, interfer- 
ometry. We stressed that the convenient environment to 
observe GPs is at zero temperature in the underdamped 
limit and that these decoherence times should be taken 
into account when extending the study of GPs to the 
nonunitary regime, especially pertinent to the applica- 
tion to fault-tolerant quantum computation [T7| . 
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